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Manipulating the energy flow of light is at the heart of modern information and communication 
technologies. Because photons are uncharged, it is still difficult to effectively control them by electrical 
means. Here, we propose a graphene plasmonic (GP) lens to efficiently manipulate energy flow by 
elaborately designing the thickness of the dielectric spacer beneath the graphene sheet. Different from 
traditional metal-based lenses, the proposed graphene plasmonic lens possesses the advantages of tunability 
and excellent confinement of surface plasmons. It is found that the proposed lens can be utilized to focus and 
coUimate the GP waves propagating along the graphene sheet. Particularly, the lens is dispersionless over a 
wide frequency range and the performance of lens can be flexibly tuned by adjusting the bias voltage. As an 
application of such a lens, the image transfer of two point sources with a separation of Ao/30 is demonstrated. 

Surface plasmon polaritons (SPPs) are electromagnetic surface waves resulting from the coupling of the 
electromagnetic field to the collective plasma excitations of metaP. Due to the ability to overcome the 
diffraction limit of light in microchip -sized devices, SPPs are considered as one of the most promising 
candidates for nanophotonic components''^. In recent years, varieties of plasmonic devices with special function- 
alities are proposed and investigated^"^. In these researches, the noble metals, such as silver and gold, have been 
popular materials to support SPP waves in the infrared and visible regimes. However, the noble metals possess 
enormous loss in the frequency regimes of interest. Especially, it is difficult to vary the permittivity functions of 
the noble metals^'^^. These insurmountable obstacles severely limit the propagation length of SPP waves and the 
tunability of the metallic plasmonic devices. 

Graphene, which is a single sheet of carbon atoms in a hexagonal lattice, has drawn much attention for its 
exceptional electrical and optical properties^ \ The complex conductivity (cTg = cTg^r + i^gd of graphene is deter- 
mined by Kubo's formula, which expresses the dependence on the radian frequency co, charged particle scattering 
rate F, temperature T, and chemical potential ^c^^'^^- Comparing with noble metals, the most advantage of 
graphene for plasmonic applications is the capability to tune surface conductivity by gate voltage, electric field, 
magnetic field, and/or chemical doping^^. In addition, the effective mode index for the TM-polarized GP waves 
along the graphene sheet can reach to ~ 100 in the mid-infrared regime^^. As a result, the GP waves can be tightly 
confined to the graphene sheet, with effective wavelength (>^gp) much smaller than free space wavelength {Xq). 
Meanwhile, the imaginary part of the wavenumber is relatively small, which means that the damping loss of GP 
waves is relatively low and the propagation length could be dozens of Xgp^^- Due to the above advantages, 
graphene has been proved to be a suitable alternative material to noble metals. Recently, a lot of GP devices have 
been proposed and investigated^^"'^. For example, Engheta et al reported that by designing and manipulating 
spatially inhomogeneous conductivity patterns across a flake of graphene, one can treat graphene as a one-atom- 
thick platform for infrared metamaterials^^. The inhomogeneous conductivity can be obtained by applying 
different gate voltages or by designing the thickness of dielectric spacer underneath the graphene sheet. Wang 
et al studied the optical coupling of GPs between graphene sheets. They employed the coupling of GPs to design 
the optical splitters, spatial switches, and Mach-Zehnder interferometers^^. Perruisseau-Carrier et al. proposed a 
terahertz plasmonic antenna by using a graphene stack. The GP antenna possesses the advantages of high input 
impedance and easy fabrication'^. All these works have proved that graphene possesses excellent confinement of 
SPPs and high tunability, which makes it an ideal material to realize various plasmonic devices working from the 
terahertz to the visible regimes 

The ability to manipulate optical fields and the energy flow of light is at the heart of modern information and 
communication technologies, as well as quantum information processing schemes'^. However, because photons 
are uncharged, it is still elusive to effectively control them by electrical means'^. In recent years, many metal-based 
plasmonic devices for subwavelength light management have been proposed theoretically and demonstrated 
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experimentally^^"^^. For example, Liu et al. numerically demonstrated 
the efficiently manipulating SPPs based on transformation optics in 
visible regime^^. They showed that the propagation of SPPs could be 
flexibly controlled by carefully tuning the dielectric material prop- 
erties adjacent to a metal. Zentgraf et al. experimentally demon- 
strated plasmonic Luneburg lens and Eaton lens by tailoring the 
topology of a dielectric layer adjacent to the metal surface^^. They 
proved that the optical properties in these lenses are changed gradu- 
ally and thus the scattering could be significantly reduced. Different 
from the above metal-based plasmonic lenses, we propose a gra- 
phene-based self-focus (Selfoc) lens for manipulating the energy flow 
of GP waves in infrared (IR) regime. By applying external gate volt- 
age, the proposed Selfoc lens can be used to focus and collimate GP 
waves that propagate along the graphene. In addition, we investigate 
the dispersion of the proposed Selfoc lens. It is found that such a lens 
is dispersionless over a wide frequency range. Finally, to illustrate the 
application of the proposed lens, we demonstrate that the image 
transfer of two point sources separated by a distance of /lo/30 (/Iq is 
the incident wavelength in vacuum) can be realized on the graphene. 
The proposed graphene-based Selfoc lens can be flexibly tuned by 
controlling the external bias voltage and possesses high confinement 



of SPPs, which paves the way for effectively manipulating energy flow 
of light at nanoscale. 

Results 

Structure model and analytical theory. As shown in Fig. 1, the pro- 
posed Selfoc lens consists of a monolayer graphene with chemical 
potential [i^ on top of a heavily doped silicon (Si) substrate separated 
by a dielectric spacer. The dielectric can be set as silicon dioxide with 
the relative permittivity of £r2 — 3.9^^'^^, and the thickness of the 
spacer h changes gradually from the position of z = 0 to periphery 
along the ± z-axis directions. The surrounding medium is assumed 
to be air and GP waves are excited and propagate along the +x-axis 
direction. A bias voltage is applied between the graphene sheet and 
doped Si substrate to change the doping of graphene by the electric- 
field effect^^'^l 

We start our study by analyzing the dispersion relation for GPs 
propagating along a monolayer graphene. In our analysis, graphene 
is treated as an ultra-thin metallic slab with a thickness of A. It should 
be noted that A is not the real thickness of graphene (—0.33 nm) and 
it is introduced only for numerical computation purpose^^'^^. This 
slab can be characterized by an equivalent permittivity of £g,eq which 
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Figure 2 | Evolution of effective mode indices of GP waves as a function of the dielectric spacer thickness with Vg = 20 V (a), and the bias voltage with h 
= 200 nm (b) for different frequencies. In calculations, the temperature Tis 300 K and relaxation time of graphene sheet is 0.5 ps. 
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Figure 3 | Lateral profile of the effective mode index Uqp (green line) and 
the thickness of the dielectric spacer h (blue line). In calculations, 
the bias voltage is 20 V and the incident frequency is 40 THz. 



is given by^^ 



:l + /o-g^o/(/coA), 



where represents the conductivity of graphene, rjo ^ 377 Qis the 
impedance of air, and ko = In/Xo is the wavenumber in vacuum. 
With Maxwell's equation, the dispersion of TM-polarized guided 
waves in graphene can be expressed as^^ 



£r2 



+ ^=0, (2) 



where /cqp is the wavenumber of the guided modes, Sj-i and are the 
dielectric permittivities of air and spacer, respectively. 

Within the random-phase approximation, the dynamic optical res- 
ponse of graphene can be derived from the Kubo's formula consisting 



of the interband and intraband contributions: = Winter + ^intm^- In 
the THz and infrared ranges, the intraband transition of electrons 
dominates. On condition that ^ ^bT, where k^ is the 
Boltzmann's constant, the surface conductivity of graphene can be 
simpUfied to a Drude-like form^^'^^ 



nh {w-\-iT~ 



(3) 



here t is the momentum relaxation time, e and H are the electron 
charge and reduced Planck's constant, respectively. The carrier relaxa- 
tion time T determined by the carrier mobility /n in graphene as t = 
fifiJ(eViy^. Recently, it has been reported that the carrier mobility of 
fi = 8000 cm^/(V-s) of graphene can be obtained by mechanical 
cleavage of bulk graphite and then transferred to Si02/Si substrate^^, 
and fi = 230000 cm^/(V-s) can be experimentally achieved in high- 
quality suspended graphene^^. When the chemical potential is 
0.15 eV with a certain gate voltage^^, the above carrier motilities corre- 
spond to the relaxation times of t = 0.12 ps and t = 3.45 ps, respect- 
ively. Here, we choose t = 0.5 ps. This value follows the ballistic 
transport features of graphene, whose mean free path was measured 
to be up to 500 nm at room temperature^^. It should be noted that the 
carrier relaxation time t can affect the imaginary part of the effective 
mode index of the graphene, which may result in the distortion of 
images and the deteriorated focusing effect. In practice, therefore, the 
carrier mobility and relaxation time should be appropriately tuned to 
reduce the distortion of images and the deteriorated focusing effect. In 
our analysis, fi^ = Wi^TUis, where Vf = 10^ m/s is the Fermi velocity 
and is the sheet doping of graphene^^. In the proposed Selfoc lens, 
can be controlled using an external bias voltage following 



(4) 



where Qx = SriSo/^ is the gate capacitance, and Vg is the external bias 
voltage^^. Based on the above equations, the effective mode indices ^qp 
{=kQp/ko) of GP waves for different thicknesses of dielectric spacer 
and bias voltages can be obtained. As shown in Fig. 2, it is found that 
Ugp can be adjusted by tuning the h and Vg. In the following discus- 
sions, we mainly employ the changing of h to control the effective 
mode index Hqp and manipulate the propagation of GP waves in 
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Figure 4 | Simulation results of the Selfoc lens, showing the amplitude of the 7 component of electric field {Ey) of GP waves. The plane GP waves are 
excited and propagate parallel (a) and obliquely (b) to the x-axis direction. Point source is illuminated at the position of z = 0 (c) and z = 150 nm (d). In 
the calculations, the frequency of the incident light is 40 THz and lateral profile of the effective mode index of graphene is the same as that in Fig. 3. 
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Figure 5 | (a) Lateral profiles of the effective mode index Uqp as a function of the incident frequencies, (b) Slope of the lateral profiles of /7gp (d/^cp/dz) for 
different incident frequencies, (c) Derivative of ^g? with respect to incident frequency (d^Gp/d/). (d) Parameter (^fre as a function of the incident 
frequencies. In the calculations, we set Vg = 20 V, t = 0.5 ps, T = 300 K and the thickness of dielectric spacer is the same as that shown in Fig. 3. 



graphene. The influences of bias voltage Vg on the performance of the 
proposed Selfoc lens are discussed in the Supplementary Materials. 

It is well-known that a Selfoc lens, which is a graded-index (GRIN) 
lens, behaves essentially like a convex lens. It can transform a cylin- 
drical wave to a plane wave^^. In general, the refractive index is high- 
est in the optical axis and decreases with transverse distance from the 
axis. The following formula closely describes the refractive index of a 
GRIN lens^^'^^ 



n{z) = no( 1- 2^'^' 



(5) 



Here, g = 2 v 2An j wis the gradient constant and An = {hq — ni)/no 

is the relative change in refractive index, Uq and Ui are the refractive 
indices at the positions of the optical axis (i.e., z = 0) and edge (i.e., 

z = ±^w) of the graphene. When the effective mode index of the 

graphene Uqp is designed and changes following Eq. (5), the gra- 
phene can act as a Selfoc lens. As shown in Fig. 2(a), one can obtain 
varies following Eq. (5) by elaborately designing the thickness of 
the dielectric spacer. Here, the refractive indices Hq and Ui are chosen 
as 87.5 and 69, respectively. According to Eqs. (2)-(5), the lateral 



profiles of n{z) and h{z) are shown in Fig. 3. It can be seen that the 
parabolic characteristic of h makes the index profile vary continu- 
ously from the optical axis to the periphery along the transverse 
direction. So when the GP waves are excited and propagate along 
the graphene sheet, they follow a cosinoidal path along the +x-axis 
direction as shown in Fig. 1(b). The period of the cosinoidal path 
("pitch" of the lens) is given by P = 2n/g, which is only determined by 
the gradient constant ^ (see Methods Section). 

Numerical results. To demonstrate the performance of the Selfoc 
lens, we employ finite- difference time-domain (FDTD) method to 
calculate the field distributions for different kinds of incident 
sources. In the calculations, graphene is treated as an ultra-thin 
metallic film with the thickness of A = 1 nm^^'^^''*\ The minimum 
mesh size equals 0.2 nm in the FDTD calculations, which is small 
enough to ensure the convergence of the results^^. For simplicity, in 
the simulation, graphene sheet is free-standing in vacuum with the 
spatially inhomogeneous conductivity patterns. Figures 4(a) and 4(c) 
demonstrate the transformation between the plane source and point 
source. When a plane (point) source is illuminated to the graphene 
sheet along the +x-axis direction, the excited GP waves can be 
transformed to be a point (plane) source. As depicted in Fig. 4(b), 
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Figure 6 | Simulation results of the Selfoc lens, showing the amplitude of the 7 component of electric field (Ey). The incident frequencies are 35 THz 
(a)-(b) and 45 THz (c)-(d), respectively. In calculations, the lateral profile of h is the same as that in Fig. 3. 



when the plane source propagates obliquely along +x-axis direction, 
it can be focused and transformed to be a point source that locates out 
of the optical axis. When the incident source is point source at the 
position of z = 150 nm, it is evolved into an obhquely transmitted 
plane source at the exit of the lens shown in Fig. 4(d). In the 
calculations, the width of the lens is chosen as w = 800 nm and 
the gradient constant can be calculated as ^ = 0.001626. So the 
pitch of the lens is P = 3864.2 nm. Considering that the proposed 
Selfoc lens is a "P/4" lens, the focal length isd' = P/4 = 966.1 nm. In 
Fig. 4, the focal length is do = 963.2 nm, which agrees well with the 
analytical results. So we can conclude that the proposed structure can 
act as a Selfoc lens and can be employed to manipulate the GP waves 
that propagate along the graphene sheet. Moreover, this Selfoc lens 
can realize the transformation optics, focusing and collimating light. 
Especially, when the proposed GP lens is used as an optical coupler, 
the full-width at half-maximum (FWHM) of the output light can be 
flexibly tuned by changing the external bias voltage. Thus, the GP 
waves can be coupled into other subwavelength optical devices with 
different dimensions (see Supplementary Materials). It should be 
noted that in the proposed Selfoc lens, the effective mode index 
ranges from 69 to 87.5 for the incident frequency of/ = 40 THz (1 
= 7.5 jim). Thus, the wavelength of GP wave is ispp ^ '^/69 = 
108.7 nm. It can be seen that the proposed graphene-based lens 
possesses excellent SPP confinement characteristics. Considering 
that the thickness of the dielectric spacer changes from 217 nm to 
350 nm, the loss of the substrate can be neglected in the simulations 
because the excited GP wave is highly confined in the graphene. 

The effective mode index of GP waves differs for varied incident 
frequencies, which means that graphene is a quite dispersive mater- 
ial. Therefore, it is necessary to analyze the performance of the pro- 
posed Selfoc lens for different incident frequencies. In Fig. 5(a), we 
calculate the lateral profiles of the effective mode index as a function 
of the incident frequencies while maintaining the dielectric spacer 
thickness unchanged. It can be seen that the effective mode index 
increases with the increase of the incident frequency, which confirms 
the dispersive characteristic of graphene. In addition, the slope of 
lateral profiles of effective mode index (dnQp/dz) for different incid- 
ent frequencies is calculated and shown in Fig. 5(b). For higher 
incident frequencies, the slope changes severely, indicating that the 
lateral profile of ^qp becomes steep and the radius of lateral profile 



curvature decreases. In Fig. 5(c), we calculate the derivative dncp/d/ 
at different lateral positions. For a specific lateral position, the deriv- 
ative is a constant, which implies that Wqp changes linearly with the 
incident frequency. We define a new parameter Cfre ^ ^Gp(/i 2:)/ 
riQpifo, z) to describe the variation of Hq^ at different positions, where 
^Gp(/j z) is the effective mode index at the position of z for the 
incident frequency of/. Here,/) ( = 40 THz) is the incident frequency 
which determines the original lateral profile of h in Fig. 3. Figure 5(d) 
shows the parameter Cfre for different incident frequencies. It is found 
that for a specific incident frequency, the parameter Cfre keeps 
unchanged at different lateral positions, which denotes that the rela- 
tive change in effective mode index An (=(^0 — ni)/no) is also 
unchanged. We can also conclude that the lateral profile of Uqp for 
different frequencies always has the form of Eq. (5), besides Uq and Ui 
have new values. Thus, though the incident frequency differs, the 
pitch of the proposed Selfoc lens is still determined by P = 27i/^ with 

g = iVlAn j w. Due to An retains constant for different incident fre- 
quencies, the pitches P and the focal length do (=P/4) are constant 
and the Selfoc lens is dispersionless. When the proposed Selfoc lens is 
designed and the dielectric spacer thickness h is fixed, the Selfoc lens 
can always focus light at the same position in the graphene sheet. This 
is because the increase of the incident frequency has two effects on 
the Selfoc lens: the increase of Uqp (as shown in Fig. 5(a)) and the 
steeper lateral profile of Hqp (as shown in Fig. 5(b)). The increase of 
Hqp results in the increase of the Selfoc lens pitch (see the Supple- 
mentary Materials), meanwhile the steeper lateral profile of Uqp 
decreases the Selfoc lens pitch^^. When these two effects are comple- 
mentary, the pitch of the proposed Selfoc lens keeps unchanged and 
the GP lens is dispersionless. As a result of this important feature, the 
GP lens designed for/o = 40 THz also works perfect well for other 
incident frequencies. 

To validate the above analyses, the field distributions are per- 
formed for different incident frequencies by using FDTD method. 
As shown in Figs. 6 (a) and 6(c), when the frequencies of the incident 
point source are chosen as 35 THz and 45 THz, the excited GP waves 
can be transformed to be plane sources. Moreover, when plane 
sources with the frequencies of 35 THz and 45 THz are illuminated, 
the GP waves can be focused and transformed to be point sources as 
presented in Figs. 6(b) and 6(d). It can be concluded that even though 
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Figure 7 | Image transfer in the graphene-based Selfoc lens, (a) Lateral profile of the input intensity of GP waves, (b) / component of electric field (Ey) of 
GP waves: the image of two point sources separated by a distance of ^o/30 ( = 250 nm) is transferred through the Selfoc lens, (c) Lateral profile of the 
output intensity of GP waves. In calculations, the incident frequency of the two point sources is 40 THz. 



the Selfoc lens is designed for 40 THz, it can still work well for other 
incident frequencies. In addition, the focal lengths for the incident 
frequencies of 35 THz and 45 THz are di = 963.1 nm and d2 = 
963.5 nm, respectively. While for the frequency of 40 THz, the focal 
length is do = 963.2 nm. It is found that the focal lengths for these 
three incident frequencies are nearly unchanged, which confirms 
that the proposed Selfoc lens is dispersionless. This important char- 
acteristic ensures that the Selfoc lens possesses a wide operating 
frequency regime, which is necessary in practical applications. It 
should be noted that the pitch of graphene-based Selfoc lens is also 
a constant for different external voltages, which is discussed in the 
Supplementary Materials. 

As an important application, we demonstrate the image transfer in 
the graphene-based Selfoc lens. As depicted in Fig. 7 (b), two incident 
point sources are separated by a distance of Xo/30 ( =250 nm) and the 
excited guided GP waves can propagate along the graphene. The 
intensity of GP waves at the output interface of the image transfer 
system is presented in Fig. 7(c). Comparing with the input intensity 
of GP waves in Fig. 7(a), the image of the two point sources can be 
well resolved at the output interface with a separation of 250.4 nm. It 
is found that the separation between the output intensity peaks only 
has 0.16% change. Especially, the separation between input intensity 
peaks is only which implies that such graphene-based lens can 
overcome the diffraction limit and realize image transfer on a mono- 
layer graphene without distortion. 

Comparing with the proposed graphene-based Selfoc lens, the 
employment of metal to realize the similar function is hardly possible. 
To get the same effective mode index as the graphene, the dimension 
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Figure 8 | Path of light propagating in GRIN media when the incident 
light is plane source (a) and point source (b). 



of the metallic devices has to be several micrometers^^ The proposed 
graphene plasmonic lens possesses ultra-compact footprint of 1.05 X 
0.85 ^m^, which is more suitable for ultra-compact optoelectronic 
devices. In addition, the characteristic of the output GP waves can 
be easily tuned by applying an external bias voltage (see the 
Supplementary Material). In contrast, it is difficult to tailor the char- 
acteristics of metallic plasmonic devices. Thus, such a graphene-based 
Selfoc lens is more suitable to manipulate energy flow of light at 
nanoscale and resolve sub wavelength features. 

Discussion 

We have proposed and investigated the GP lens for effectively 
manipulating energy flow. The effective mode index of GP waves 
can be easily tuned by changing the dielectric spacer thickness or the 
external bias voltage. Based on this principle, the graphene-based 
Selfoc lens is elaborately designed and FDTD simulation is per- 
formed to characterize the optical properties of the Selfoc lens. It is 
found that the proposed Selfoc lens can be used to collimate and 
concentrate GP waves at subwavelength scale. We also demonstrate 
that the proposed Selfoc lens is dispersionless over a wide frequency 
range. As an important application of the graphene-based Selfoc lens, 
we demonstrate that the image of two point sources separated by a 
distance of XqITd can be well resolved on the graphene. Comparing 
with metallic plasmonic devices, such a graphene-based Selfoc lens 
possesses ultra- compact dimension and flexible tunability, and 
hence is more suitable for manipulating energy flow at nanoscale 
and resolving subwavelength features. 

Methods 

Light propagation in GRIN media. In the proposed structure, the different substrate 
thicknesses result in the spatially inhomogenous effective mode index across a flake of 
the graphene. The index profile varies continuously from the optical axis to the 
periphery along the transverse direction. In our studies, the mode index is designed 
and varies quadratically with the distance in the transverse direction and the graphene 
is regarded as Selfoc lens^'''^^. 

According to the Fermat's principle, the light equation for the light propagating in 
GRIN media can be expressed as^'^. 



dT 



(6) 



where ds is a differential arc length along the ray trajectory, is the position vector 
and W nis the gradient of n. In the xz-plane, with considering the paraxial light, ds ~ 
dx and Eq. (6) is: 



d^ 



dn 
Tz' 



§n dz 
dxdx' 



(7) 



The refractive index of the medium varies in the z-axis direction satisfies Eq. (5) and - 
g^z^ <C 1, thus Eq. (7) can be simplified as: 
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the general solution of Eq. (8) is 

z = B cos{gx) + C sin(^x) , 



(8) 



(9) 



where B and C are constants determined by the initial conditions. According to Eq. 
(9), the slope of light line can be derived as, 



dz 



- Eg sin (gx) + Cg cos (gx) . 



(10) 



Here, the initial conditions are assumed as (Kq, Zq), which means that the slope and 
position of the incident light is Kq and Zq, respectively. By solving Eqs. (9) and (10), the 
constants B and C are calculated as 5 = Zq, C = Ko/g, then Eq. (9) can be obtained as. 



= Zo cos{gx) + (Ko/g) sm{gx). 



(11) 



When the incident light at the position of Zq is plane wave parallel to the optical axis 
(i.e., the initial conditions are Kq = 0 and z = Zq), according to Eq. ( 1 1), we can get the 
light equation 



z = zocos(gx). 



(12) 



which means that the light propagating in GRIN media follows a cosinoidal path and 
the period (also called "pitch") of the Selfoc lens is P = 2n/g. 

As the examples, the path of the light propagating in GRIN media is shown in Fig. 8. 
It can be seen that when the proposed Selfoc lens is used as an optical coupler, the 
length of the lens should be P/4. When it is used to realize the image transfer, the 
length of the lens should be P/2. 

When the effective mode index has an increment of An, by solving the above 
equations, the pitch of the Selfoc lens can be easily obtained as. 



^, In ln-\-An 



(13) 



It can be seen that the pitch of the Selfoc lens increases with the increase of effective 
mode index. 
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